Optimal estimation of squeezing 
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We present the optimal estimation of an unknown squeezing transformation of the radiation field. 
The optimal estimation is unbiased and is obtained by properly considering the degeneracy of the 
squeezing operator. For coherent input states, the r.m.s. of the estimation scales as {2\/n)~ 1 versus 
the average photon number n, while it can be enhanced to (2n) _1 by using displaced squeezed 
states. 

PACS numbers: 



Squeezed states are characterized by a phase- 
dependent redistribution of quantum fluctuations such 
that the dispersion in one of the two quadrature com- 
ponents of the field is reduced below the level set by the 
symmetric distribution of the vacuum state or a coherent 
state . Such a property has been used to raise the sen- 
sitivity beyond the standard quantum limit [2, |3| an d to 
enhance interferometric [4| and absorption measurements 
H , along with optical imaging applications 0, @ . 

Even though squeezed states have been studied exten- 
sively during the last three decades, the attention to the 
problem of estimating an unknown squeezing transforma- 
tion is relatively recent, and very few results are known 
about it. The first attempt to quantify the accuracy lim- 
its imposed by quantum mechanics was presented in Ref. 
0, in the case of squeezing in a fixed direction. Here, the 
squeezing transformations form a one-parameter group, 
and the estimation problem is closely similar to the prob- 
lem of phase estimation 0, (for this reason the name 
hyperbolic phase estimation has been also used). The 
basic idea underlying the estimation strategy is to find 
a measurement that projects the quantum state on the 
vectors that are canonically conjugated via Fourier trans- 
form to the eigenstates of the squeezing generator. How- 
ever, as we will show in this Letter, the scheme of Ref. 
is not optimal and is biased: neither the mean value 
nor the most likely one in the probability distribution co- 
incide with the true value of the squeezing parameter. In 
other words, the estimation is biased, and the presence 
of such a bias suggests that the proposed scheme is not 
optimal. 

More recently, the estimation of squeezing has been 
considered in connection with cloning In this case 
the unknown squeezing is estimated from a number of 
identical copies of the same unknown squeezed state. 
However, this approach does not work when only a sin- 
gle copy is available, and the problem of the bias and the 
optimality of the estimation is left open. 



In this Letter, we will present the optimal estimation 
of an unknown squeezing transformation in a given direc- 
tion, acting on an arbitrary state of the radiation field. 
This problem is the optimal estimation of squeezing in 
an experimental situation where a degenerate paramet- 
ric amplifier is pumped by a strong coherent field with a 
fixed phase relation with the state to be amplified. 

We will show that the optimal measurement is unbi- 
ased, provided that one properly takes into account the 
degeneracy of the squeezing operator. Due to such a de- 
generacy, the Fourier transform of the eigenstates of the 
squeezing operator is not uniquely defined, and, in order 
to obtain the best estimation strategy, one has to perform 
an optimization similar to that of phase estimation with 
degeneracy 0. Accordingly, the optimal estimation of 
squeezing depends on the chosen initial state of the ra- 
diation field. Also, the optimization performed here is 
analogous to that of Ref. in the case of estimation 
of rotations, namely it properly takes into account the 
equivalent representations of the group of parameters. 
In fact, the degeneracy of the squeezing operator corre- 
sponds to the presence of equivalent representations of 
the related one-parameter group. 

We will derive our results in the framework of quantum 
estimation theory 0, 0] , upon defining optimality as 
the minimization of the expected value of a given cost 
function, which quantifies the deviation of the estimated 
parameter from the true one. According to the minimax 
approach, the optimal estimation strategy will be the one 
that minimizes the maximum of the expected cost over all 
possible true values of the unknown squeezing parameter. 
In analogy with the class of cost functions introduced 
by Holevo 0, 0] for the problem of phase estimation, 
we introduce here a class of cost functions including a 
large number of optimality criteria, such as maximum 
likelihood, and maximum fidelity. We will show that our 
estimation strategy is optimal according to any function 
in such a class. 
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In the following, we consider a single-mode radiation 
field with bosonic operators a and a', satisfying the 
canonical commutation relations [a, a^} = 1. The squeez- 
ing operator is defined as follows 



S(r) 



exp 



2 [a 



f2 



(1) 



where r is a real parameter. Given a pure state of 
the radiation field, we want to find the optimal mea- 
surement that allows one to estimate the parameter r 
in the transformation \%b) — ► S , (r)|-0). In the quadra- 
ture representation ^>{x) — (x\th), where |x) denotes the 
Dirac-normalized eigenvector of the quadrature operator 
X = (a + a+)/2, the effect of squeezing on the wavefunc- 



tion is given by ip(x) 



-r/2 



The squeezing operator can be written as S(r) 
, where K is the Hermitian operator K = i(a 



-irK 



f2 



a )/2, that generates the one-parameter group of squeez- 
ing transformations. The spectrum of the generator K is 
the whole real line, and the eigenvalue equation reads 



K\fi,s) = (J,\fi,s) 



(2) 



where /i 6 R is the eigenvalue, and s is a degeneracy in- 
dex with two possible values ±1. The explicit expression 
of the generalized eigenvectors of K in the quadrature 
representation is given by 0] 



1 



(x\n, s) = 



'2tt 



(sx) , 



(3) 



where 9(x) is the Heaviside step-function [8(x) = 1 for 
x > 0, 9(x) = for x < 0]. The vectors s) are orthog- 
onal in the Dirac sense, namely (fi, r\v, s) = S rs S(fi — is), 
and provide the resolution of the identity 



dfj, IL. = 1 , 



(4) 



where 11^ = X} s =±i Ia 4 ' s )(^' s I ^ s ^ ne projector onto the 
eigenspace of K corresponding to the eigenvalue /i. 

Let us denote by 7i M the two-dimensional vector space 
spanned by \ fi, ±1). In this complex vector space, we can 
consider the usual scalar product and the corresponding 



norm, namely if = J2 S = 
r K lll then its norm is 



±i 



s) is an element of 



= (£« 



^| 2 ) 1/2 . Using 



the completeness relation (0J , we can write any pure state 
lib) eH as 



where c p = ||ILJ 



|, and 



|n M |V)ll 



(5) 



(6) 



is the normalized projection of \ib) onto H^. The rep- 
resentation of a state as in Eq. corresponds to the 



fact that the Hilbert space Ti. can be decomposed as a 
direct integral Tt = f^^dfj. Ti^. In this representation 
the effect of a squeezing transformation is given by 



S(rM) = 



+oo 



dfi Cfj 



-IT fX 



(7) 



i.e. the squeezing operator introduces a different phase 
shift in any space TC^. Notice that the states (0) all lie 
in the subspace 



dfj, v^ib^) I Vf, e L 2 (R) 



(8) 



The problem of squeezing estimation in the representa- 
tion JJJ becomes formally equivalent to the problem of 
phase estimation. 

In order to optimize the estimation of squeezing, we de- 
scribe the estimation procedure with a positive operator 
valued measure (POVM) P(f). The probability distribu- 
tion of estimating f when the true value of squeezing is r 
is then given by p(f\r) = Ti[P(f)S r pSl]. The optimality 
criterion is specified in terms of a cost function c(f — r), 
that quantifies the cost of estimating f when the true 
value is r. Once a cost function has been fixed, the opti- 
mal measurement is defined in the minimax approach as 
the one that minimizes the quantity 



max 

rGK 



df p(f\r)c(f — r) 



(9) 



namely the maximum of the expected cost over all pos- 
sible true values. Generalizing the class of cost functions 
introduced by Holcvo for phase estimation 0, we con- 
sider here cost functions of the form 



c(r) 



+ OC 



dfi a M cos(^tr) 



(10) 



where a M < for any /i > 0. This class contains a 
large number of optimality criteria, such as the maxi- 
mum likelihood Cml{t) = —8{r), and the maximum fi- 
delity c F (r) = l-\(iJj\S(r)\i;}\ 2 . 

Due to the group symmetry of the problem, instead of 
searching among all possible measurements for optimiza- 
tion, one can restrict attention to the class of covariant 
measurements which are described by POVMs of 
the form P(f) = S{f)£,S{f)\ with f > such that 



dr S(r)£S(r)* 



1 



(11) 



The probability distribution p{f\r) related to a covariant 
measurement will depend only on the difference f — r, 
and this means that the estimation is equall y g ood for 
any possible value of the unknown squeezing 15]. 

The optimization of the covariant measurement for any 
cost function in the class (|10fl can be obtained as in the 
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case of phase estimation with degeneracy 9]. The opti- 
mal covariant POVM is then given by 



P(r) = \v(r))(v(r)\ 



where 



\r)(r)) 



+00 



v/27T 



(12) 



(13) 



Notice the correspondence of \r](r)) with the vectors 
\e(tp)) = X^^Lo ^/2^ \ n ) tnat ar i se m the context of op- 
timal phase estimation (here \n) are the non-degenerate 
eigenvectors of the photon number operator at a). The 
vectors \f](r)} are orthogonal in the Dirac sense, namely 
the optimal POVM is a von Neumann measurement. The 
projection j^) of Eq. JBJl in the expression of \r](r)} 
makes the optimal measurement depend on the input 
state IV'). Accordingly, one obtains non-commuting ob- 
servables, corresponding to different input states. The 
normalization of the POVM l|12|) can be easily checked, 
since /_ °° dr P(r) = 1^ , where 1^, is the identity in the 
subspace Tt^ defined in Eq. (jHJ. Clearly, the P(r) can be 
arbitrarily completed to the whole Hilbert space, without 
affecting the probability distribution of the outcomes. 

Using Eq. 10, the optimal probability distribution for 
an input state \ip) is given by 



p(f\r) = (ip\S(r^ P(f) S(r)\ip) 



1 

2^ 



-i(r—r)fi 



.(14) 



Since the probability distribution depends only on the 
difference f — r, from now on we will write p(f — r) 
instead of p(r\r). Representing the projection 11^ as 
IIu = l^e^O--^, the probability distribution of 

Eq. (|14|l can be rewritten as 

'■+ 00 dfi 



p{r) 



+ OO 



-irfi 



2vr 



d\e iX » (ip\S(\)\ip) 



(.15) 



The optimal measurement l|12[) can be compared with 
that given in Ref. |jj , which is described in our notation 
by the POVM 



Hr)= E \Vs(r))(Vs(r)\ 



s=±l 



where 



Mr)) = 



+00 



dfj, 
J2^ 



(16) 



(17) 



Using Eq. it is easy to see that \r]±(r)) are eigen- 
vectors of the quadrature X corresponding to the eigen- 
values ±e r , and hence the POVM (|16|) corresponds to 
measuring the observable In \X\, independently of the in- 
put state. The measurement P(r) is not optimal, and 



gives a biased probability distribution, namely the aver- 
age value of the estimated parameter does not coincide 
with the true value, and also the most likely value in the 
probability distribution is not the true one (see, e.g., the 
asymmetric probability distribution for the vacuum state 
in Fig. 1). Such drawbacks do not occur in the optimal 
probability distribution (|14|) . Notice also that the mea- 
surement P(r) is "rank-two" in the subspace TL^p of inter- 
est, while the optimal measurement is "rank-one". The 
differences between the two measurements can be under- 
stood intuitively as follows. Essentially, both POVMs are 
based on the Fourier transform of the eigenvectors of the 
operator K. However, since the Fourier transform is not 
uniquely defined due to the degeneracy of K , one should 
optimize it versus the input state. 

p (r) 




FIG. 1: Probability distributions for the estimation of squeez- 
ing on a vacuum input state. The asymmetric distribution 
comes from the suboptimal measurement of Ref. in Eq. 
11611 . The symmetric distribution corresponds to the optimal 
measurement of Eq. W'll . 



In the case of a coherent input state \a), the probability 
distribution (|15|) can be specified as follows 



p(r) 



+30 



dfj, 
2tt" 



(18) 



i> + 00 



dX 



V cosh A 



I Q I 2 

= iA/i gj tanh A(a* 2 — a 2 ) g ^L-l ^ 



This probability distribution has been plotted for increas- 
ing real values of a in Fig. [21 where one can easily observe 
the corresponding improvement in the estimation. 

For large values of \a\, from Eq. I|19f> one obtains 
asymptotically the Gaussian distribution 



p{r) 



2d 



(19) 



that provides a r.m.s error on the estimation of r as 
Ar = 1/(2%/^), where n = |a| 2 is the mean photon 
number. This scaling can be improved to Ar = l/(2n) 
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P (r) 




FIG. 2: Optimal probability distribution of squeezing for in- 
put coherent states. The distribution becomes sharper for in- 
creasing values of the coherent-state amplitude (a = 1, 2, 4.) 



by using displaced squeezed states \a, z) — D(a)S(z)\0), 
with a, z £ R. In fact, from the relation D(a)S(z) = 
S(z)D(ae~ z ), the probability distribution p(r) is given 
by Eq. I|18|) just by replacing a with ae~ z . In the asymp- 
totic limit of large number of photons n — \a\ 2 + sinh 2 z, 
the minimization of the r.m.s. gives the optimal scaling 
Ar = l/(2n), for a = y/nj2 and z = -l/21n(2n), and 
this corresponds to approximate the eigenvectors of the 
quadrature operator X. 

In the asymptotic regime, the optimal performance can 
be achieved simply by measuring the quadrature X and 
estimating f = ln|x/o;|, in correspondence to the out- 
come x. However, it is important to stress that homo- 
dyne measurement becomes optimal only for particular 
input states and in the asymptotic limit of large energy, 
while for finite energy the optimal measurement is de- 
scribed by the POVM in Eq. lfT2l . 

In conclusion, we presented the covariant measurement 
for estimating the squeezing that is optimal for a large 
class of figure of merit. The optimal detection is given 
by a suitable Fourier transform of the eigenstates of the 
generator of squeezing. In fact, due to the degeneracy 
of the squeezing operator, there is a freedom in choosing 
how to perform the Fourier transform, and the choice 
must be optimized according to the input state. Hence, 
for different input states one has different optimal esti- 
mations corresponding to different observables. The op- 
timal measurement leads to an unbiased estimation, and 
the outcome of the measurement that is most likely to be 
obtained coincides with the true value of the unknown 
squeezing. For coherent input states the r.m.s. error 
scales as l/(2\/E) with the number of photons, while for 



displaced squeezed states one achieves 1 / (2n) scaling. In 
the asymptotic regime, such a scaling can be obtained ex- 
perimentally by homodyne measurement. The presented 
scheme applies to the problem of optimal characteriza- 
tion of nondegenerate parametric amplifiers. 
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